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Abstract 

Yang-Mills gauge theory models on a cylinder coupled to external matter 
charges provide powerful means to find and solve certain non-linear integrable 
systems. We show that, depending on the choice of gauge group and matter 
charges, such a Yang-Mills model is equivalent to trigonometric Calogero-Moser 
systems and certain known spin generalizations thereof. Choosing a more gen- 
eral ansatz for the matter charges allows us to obtain and solve novel integrable 
systems. The key property we use to prove integrability and to solve these sys- 
tems is gauge invariance of the corresponding Yang-Mills model. 



1 Introduction 

In a previous paper [ij] we presented a novel class of integrable spin-particle systems 
which contains known integrable systems of Calogero-Moser (CM) type and certain 
known spin generalization p, ^ thereof as special cases, and many other generaliza- 
tions which (to our knowledge) were not known before. Our method not only allowed 
us to find and prove integrability of these models but also to solve them explicitly. 
Subsequently two alternative derivations of these models also proving integrability 
were given by Polychronakos 0, 0. 

In the present paper we give a more detailed account of our approach and also 
present extensions of our previous results. We made some effort to make the paper 
easily accessible to different readers: those who quickly want to get the flavor of how 
our method work, but also those who are interested in the details. In the rest of this 
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section we give some introduction to CM type systems, describe our method, and 
then give a plan for the paper. 

About thirty years ago it was discovered that a dynamical system of particles on 
the line described by the Hamiltonian 

is completely integrable when the pair potential v{r) equals a Weierstrass elliptic func- 
tion p(r), important special cases of which are a? / sin^ (ar), and a? / sinh^ (ar) 
[0] (for an early review see |^). Subsequently these models have received much in- 
terest in different contexts, and various generalizations of these models have been 
found and studied (for recent reviews see, e.g., |§, |^). Recently it was observed that 
these models can be obtained from Yang-Mills gauge theories on the cylinder coupled 
to particular non-dynamical matter charges Exploring this relation further we 



found and solved a large class of novel integrable spin-particle systems In this 
paper we give a more detailed account of these and some new results. The key prop- 
erty we use to prove integrability and to solve these systems is gauge invariance of 
the corresponding Yang-Mills model. This allows us to use different gauges, i.e., to 
impose different constraints compatible with gauge invariance. There is a gauge in 
which the dynamics of such a Yang-Mills model is equivalent to the dynamics of a 
CM type system, whereas in another gauge the dynamics is free and the solution can 
be found trivially. Thus the solution of the former system can be obtained from the 
latter solution by constructing a certain gauge transformation. 

Relating a CM type system to a Yang-Mills theory makes integrability obvious 
and the construction of conservation laws (nearly) trivial. Moreover, this relation 
provides a rather simple method to construct an explicit solution. (This method can 
be regarded as an extension of the projection method 0.) Since there is a large 
freedom in choosing the external matter charges, one can obtain and solve a large 
number of different integrable systems. We believe that it should be possible to find 
other integrable models using our method, for example by considering more general 
gauge groups etc. 

The plan of this paper is as follows. In Section ^ we summarize the facts about 
Yang-Mills gauge theories on a cylinder which we need in the sequel. The derivation 
of a certain class of dynamical systems from these gauge theories is explained in 



^This is actually implicit already in earlier work; see, e.g., [O 
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Section ^. We present two different arguments: the first argument (Section is 



quick but only heuristic, and the second (Section |3^) is somewhat less intuitive but 
rigorous. In Section ^ we show how to exploit the relation of these dynamical systems 
to gauge theories to exphcitly solve interesting special cases of these systems, and 
we also show how to obtain Lax pairs and conservation laws in our formalism. As a 
warm-up, we first show how known results about CM models with potentials 
a} j sin^(ar) and c? j sinh^(ar) can be recovered (Sections |4.1| -4.5). We will also refer 
to the former Calogero model and the latter as Sutherland model. We then 

derive and extend the results for the novel systems found in |l| (Section ^!4|) . A novel 
class of models which can not be solved in such an explicit manner but still should 



be integrable is discussed in Section ^ 



Readers who only want to get a flavor of how our methods work are advised to 



read the beginning of Section H and proceed to Sections 4.2 and 4.4 



Notation. We denote as gl^y the complex N x N matrices, GLjv the complex invert- 
ible N X N matrices, and / the N x N unit matrix. 



2 Yang-Mills theory on a cylinder 

In this section we fix notation and define the gauge theory models of interest for our 
purposes. We also discuss different gauge conditions which will play important roles 
in the sequel. 

2.1 Definition 

We consider Yang-Mills theory on a cylinder with external matter sources. We restrict 
ourselves to a two-dimensional spacetime which is a cylinder, i.e., the time coordinate 
is t G M, and space is a circle parameterized by x G [— 7r,7r]. In the following, 
/i, G {0,1} are spacetime indices, and repeated spacetime indices are summed 
over.Q Our metric tensor is diag(l, —1). 
Our starting point is the Lagrangian 

C{t, x) = -i- Qtr [F,,{t, x)F^''{t, x) + A,{t, x)f{t, x)]^ (2) 
^Note that we will use this summation convention only for spacetime indices. 
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where we set 

f^p. (3) 

This Lagrangian describes a dynamical Yang-Mills field A^, coupled to an external 
matter current j'^ whose spatial component vanishes. We shall denote the temporal 
component j° = p of this matter current as charge. The Yang-Mills curvature F^^, is 
defined as 

F^,u = d^A^ - d^A^ + ig[A^, A^\. (4) 
Note that the only non-trivial component of this is 

E{t,x) ■= Foi{t,x) = -Fio(t,x). (5) 

We assume a representation of the structure groups such that A^ and p are gljy- valued 
functions. 

2.2 Gauges 

The Lagrangian defined in Eqs. and is obviously invariant under all transfor- 
mations 



Af,{t,x) A^{t,x) = U~^{t,x)Af,{t,x)U{t,x) + —U'\t,x)d^U{t,x) 
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p{t,x) p^{t,x) = U ^{t, x)p{t, x)U{t, x) 

E{t,x) ^ E^{t,x) = U-\t,x)E{t,x)U{t,x) (6) 

for all differentiable GL^r-valued functions U on spacetime (the last relation is of 
course not independent since it follows from the first). This implies that the equa- 
tions of motion derived from the Lagrangian are also invariant under all these 
transformations. We will consider functions U {t, x) which are differentiable only lo- 
cally in time, i.e., it is differentiable for t in some non-empty interval (^1,^2) and all 
X G [— TT, tt]. We denote the set of all such functions as G{ti,t2)- We refer to the trans- 
formations in Eq. @ as gauge transformations. Gauge invariance, i.e., invariance 
under all gauge transformations, will allow us to impose additional constraints on 
the gauge fields. We call these constraints gauges. There are two different gauges 
which will be important for us. 



4 



Diagonal Coulomb Gauge 

This gauge is defined by the condition that Ai is a diagonal matrix which is inde- 
pendent of X, i.e., 

Ai(t, x) = Q{t) = diag {q\t), q\t), . . . , . (7) 

To show that this indeed is a gauge one has to prove that it is always possible 
to find a gauge transformation bringing the Yang-Mills field into the form of Eq. 
(0), i.e., for each (generic) Yang-Mills configuration Aiit^x) one can find a gauge 
transformation U such that A^[t,x) defined by Eq. has the form of Q{t) in Eq. 
(1^. For completeness we now recall a proof of this by explicit construction of the 
required U (see |]T| and references therein). 

We first note that the condition A^ = Q is equivalent to the differential equation 
diU + igAiU = igUQ. To solve this equation we first consider the boundary value 
problem 

diS{t,x) + igAi{t,x)S{t,x) =0, S{t,-n) = I. (8) 

This problem has a unique solution (see, e.g., Ref. [|T^, Theorems 1.1 and 2.1) which 
we write as 

S{t,x) = Pexp (^-tg £jyA,{t,y)^ (9) 

where the symbol Pexp denotes the path ordered exponential.^ It is important to 
note that S{t, x) is not a gauge transformation since it is not periodic in x (its values 
at a; = — vr and vr are different in general). However, it can be used to construct a 
gauge transformation as follows. We first introduce an important technical condition. 
We call a gauge field Ai{t, x) regular (at time t) if the corresponding matrix S{t, n) is 
non-degenerate, i.e., all the eigenvalues of S{t,7i) are different. If Ai{t,x) is regular 
then there exists an invertible matrix V{t) diagonalizing S{t, vr) (see e.g. Theorem 
10.2.4. in 0), i.e., 

Vit)-^Sit, n)Vit) = e-^52-QW (lo) 

for some diagonal matrix Q{t) = diag(g^(t), . . . , g^(t)). We now claim that the 
function 

U{t, x) = S{t, x) V(t)e'^'("+'^)«W (11) 

is periodic. Indeed, 

U{t,7r) = V{t)V{t)-'S{t,7r)V{t)e''^''^^'^ = V{t) = U{t,-ix) 
■^in the terminology of Ref. fl^, S{t,x) is identical with the product integral H-tt e~*^^i 
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where we inserted V{t)V(t)^^ = I and used Eq. (]I0|). Moreover, U Eq. ( [Tl| ) satisfies 
diU + igAiU = igUQ equivalent to = Q. It is also easy to see that U(t,x) is 
invertible (see e.g. Ref. [|T^], Section 1.1). We are left to show that U{t,x) is a differ- 
entiable function on spacetime. This might seem trivial but is actually not since V{t) 
and Q(t) as defined in Eq. (|10D can be discontinuous in t; see, e.g., [ITSf. However, if 
Ai{t, x) is regular for all t in an open time interval (ti, ^2) then Q{t) and V{t) can be 
chosen to be differentiable in this time interval DTSf, and U{t, x) in Eq. (|TlD is indeed 
a differentiable function, i.e., a gauge transformation in G{ti,t2)- We call a Yang- Mills 
field Ai{t, x) generic if it is regular for all times t G K. 

It is worth noting that the function V{t) is not unique: all transformations 



with arbitrary differentiable and non-zero functions (^{t), a = 1,2, . . . , N, are com- 
patible with the conditions determining V(t). This corresponds to the residual gauge 
freedom which remains after imposing the diagonal Coulomb gauge: gauge trans- 
formations Eq. d^) with U{t,x) = D{t) (diagonal and independent of x) leave the 
diagonal Coulomb gauge Eq. (|^ invariant but act non-trivially on Aq and p. 

We note that the Eqs. (H) and (|10D provide the recipe how to compute Q{t) from 
a given Yang- Mills configuration Ai[t,x). This will play an important role for us in 
the sequel. 

Remark: We note that in our examples later the eigenvalues e"*^^'^^"*^*-' of S{t, tt) 
have the following physical interpretation: q°'{t), a = 1,...,N, correspond to the 
positions of interacting particles in a dynamical system. The technical condition of a 
Yang-Mills field being regular in some interval thus means that the particles do not 
collide with each other. Since the particle interactions are repulsive in all dynamical 
systems we encounter, it is plausible that regularity holds for all times, but we will 
be able to prove this only in certain special cases. In general we can prove our results 
only locally in time, i.e., for time intervals where no particle collisions occur. 

Weyl gauge 

This gauge is defined by the condition 



V{t) V'{t) = V{t)D{t), D{t) = diag{d\t), d^{t)) 



(12) 



Ao{t,x) = 0. 



(13) 
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To show that this is a gauge one can use a similar argument as above, i.e., for a given 
Ao(t,x), the function 

W{t, x)=V exp (^-ig dt' A^it' , x)^ , (14) 

is such that A^{t^x) Eq. (D vanishes. Note that W{t,x) here is indeed a gauge 
transformation (i.e., a periodic, invertible, and differentiable function on spacetime). 
We note that if our spacetime was a torus and not a cyhnder, the Weyl- and the 
diagonal Coulomb gauge would be more similar: one of them is obtained from the 
other by interchanging the space- and time coordinates x and t. The reason why the 
Weyl gauge is simpler on the cylinder is that there is no periodicity condition in t. 



3 Dynamical systems from Yang-Mills theories on 
the cylinder 

In this section we show how to derive certain dynamical systems from Yang-Mills 
systems on a cylinder with non-dynamical external matter charges. We will use two 
different arguments which both lead to the same results. In Section |3.1| we use a 
canonical procedure to formulate our Yang-Mills system as an infinite dimensional 
Hamiltonian system with a constraint called Gauss' law. By solving Gauss' law in 
the diagonal Coulomb gauge Eq. (^ we obtain an integrable non-linear Hamiltonian 
system which generalizes the Calogero- and Sutherland models discussed in Section 
m. This Hamiltonian method is conceptually simple but not quite satisfactory from 
a mathematical point of view. To make it mathematically rigorous would require a 
deeper analysis which is beyond the scope of the present paper. Instead, we discuss 
an alternative and complimentary method in Section |3.2| . This method uses only the 



equations of motion and thus circumvents all the mathematical difficulties which one 



would have to face in a rigorous discussion of the method in Section 3.1 



3.1 Hamiltonian approach 

In the following it is useful to denote the elements of the matrices M = Afj_,E,p 
by M"^, a,f3 = 1,2,...N, and tr the usual N x N matrix trace. Noting that 
tr(MM') = Y,a 13=1 ^"'^ following the standard canonical procedure [|16 



we obtain from the Lagrangian (E])-(H) (up to a surface term which vanishes due to 
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the periodicity of AqE in x) 



Ti- = — I dx 



tr[^E^-Ao[d,E + ig[A,,E]-p] 



(15) 



The variable conjugate to Ai^{x) is E^"-{x), which imphes the Poisson brackets 

{Af (x), E^'^\y)] = 27r5"'3'5^"'5(x - y). (16) 

Since the Lagrangian is independent of Aq^[x) its conjugate momentum, which we 
denote as no°'(x), has to vanish: Ho — 0. Here '~' means that this is to be regarded 
as a constraint [0. This primary constraint imphes the secondary constraint OqU-q = 
{no,i/} ~0, i.e., 



G{x) := diE{x) + ig[Ai{x), E{x)] - p{x) ^ 0, 



:i7) 



which is called Gauss' law. The tertiary constraint doG{x) = {G{x), H} ^ together 
with {Ai{x),p{y)} = {E{x), p{y)} = then fixes the Poisson bracket of the charges 
(see Appendix ^ for details) 



{p"^(x),p"'^'(l/)} ~ ig27r [p°^'(x)5^"' - p"''^(x)(5'^'"l S{x - y). 



:i8) 



Note that, 



{G°^(x),G'°'^'(y)} ~ tg27T [G"^'(a;)5^"' - G°''^(x)5^'"] 6{x - y), 

and the Poisson brackets above also fix the time evolution of the charges, dop{x) = 
{p{x),H}. 

We now exploit the gauge freedom and impose the gauges discussed in the last 
section. We shall do that by replacing '~' by '=' in the equations given above, i.e., 
we ignore possible subtleties and treat constraints like strict equalities. We therefore 
regard the following as a simple but only heuristic argument, as discussed above. 



Diagonal Coulomb Gauge 

As shown in the last section, it is always possible to transform Ai to a x-independent 
diagonal matrix. This means that the Yang-Mills field only has a finite number of 
true dynamical degrees of freedom. In order to obtain a Hamiltonian for these degrees 
of freedom only, we impose the gauge condition Eq. (|^). We then use Gauss' law to 
eliminate all but the true dynamical degrees of freedom. 



8 



To do this we will use Fourier transformations 

E^^n) = r rfa;e-*"^'E"^(x), nez (19) 



and similarly for Ai and p. Using this and imposing Eq. the Gauss' law can be 
rewritten in component form as 

+ =p"^(n). (20) 

For a = P and n = 0, ( pOD implies 

p°°(0) = 0. (21) 

In all other cases, Eq. (|20|) can be solved for E°'^{n). Inserting these in the Hamilto- 
nian (P^, we obtain 



1 /-TT 1 „ ^ 

n = - dxtr{E{xr) = —j: E-'i^)E'''{-n) = 

N E--(^0)E--{0) 1 ^ p-P{n)pP-{-n) 
h 2(2vr)2 ^-'kii: "° \n + g[q<^-q^]r 

Here a comment on our notation is in order: for a = (3 and n = the last term seems 
ambiguous since it formally is 0/0. However, from our derivation above it is clear 
that this has to be interpreted as 0, i.e., here and in the following 

(1 - 5„,o5"'^)(- ■ •) = if n = and a = /3 

even if (■ ■ ■) happens to be infinite. 

We now find it convenient to introduce the following suggestive notation, 

^ (22) 

/TT 

Moreover, we will see later that the functions ylo"(t, 0) = a°'{t) play a special role 
and can be chosen arbitrarily. Thus despite of the constraint Eq. (^1]) we can leave 
the corresponding term oc Yl,a ci°'{t)p°'°'{t, n = 0) in the Hamiltonian. We thus obtain 

^=E^ + AEE(1- Sn,05'^') ^\^Ta Z + A E «"P""(0) (23) 
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where a"(t) play the role of an arbitrary external time dependent field which con- 
tribute to the time evolution of 

Using Eq. (|T^) it is easy to verify that 



) [ 27r 

Hence, it is natural to interpret p"" as particle momenta and as the correspond- 
ing position variables. From the Poisson bracket (|T^) we obtain the corresponding 
Poisson bracket for p, 

{p"^(n),p"'^'(m)} = z(727r[p"^'(n + m)5^"' -p"''^(n + m)5^'"]. (25) 

We interpret the p"^{n) as spin degrees of freedom. Hence, we have obtained a dy- 
namical system described by the Hamiltonian (^) , together with the Poisson brackets 



( |24D and p5|). All other Poisson brackets vanish. 

It is interesting to note that the Hamiltonian Eq. ( p3|) is invariant under the 
transformations 

q^^q'^^ , p"^(n) -^p"^(n + m° -m^) (26) 

9 

for all integers nf. Thus, if g"", a = 1, . . . , A^, are real, this model describes particles 
moving on a circle of length 1/g and interacting with a potential whose strength 
depends on the spin degrees of freedom. Moreover, this Hamiltonian is invariant 
under all permutations of the particle labels. It is worth noting that the existence 
of these symmetry transformations is due to the Gribov ambiguities of the diagonal 



Coulomb gauge as discussed in Ref. [|T7 



Using the Hamiltonian Eq. (|2^) and the Poisson brackets given above it is straight- 
forward to derive the equations of motion d^X = {X, 7i} for X = q°',p'^ and p"^{n). 
(We will write down these equations in Section |3]^.) Note that these equations are 
highly non-linear. 

Weyl gauge 

In the Weyl gauge, the Hamiltonian is simply 

n = ^ r dxtr{E{xf) (27) 

477 J-TT 
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which shows that we have a free (non-interacting) system. Still, the system is not 
trivial since we have to account for Gauss' law Eq. i\rf\) . However, it is easy to see 
that all quantities tr[E{x)"'] (arbitrary x G [— vr, vr] and n G N) are conserved in time 
t. Thus there is an infinite number of conservation laws. (Note that we do not 
claim that these conservation laws are independent.) Since the system is free, we can 
actually solve the system in the Weyl gauge (this will be done in the next section). 

3.2 Alternative approach based on equations of motion 

Our discussion in the last section suggests that our gauge theory is an integrable sys- 
tem and equivalent to a non-linear dynamical system which is obtained by imposing 
the diagonal Coulomb gauge. Instead of trying to make the argument in the last 
section rigorous (which would be interesting but is beyond the scope of the present 
paper) we now present an alternative, less intuitive, argument leading to the same 
conclusions but avoiding the mathematical difficulties of the canonical procedure. 

We start with the Lagrangian equation obtained from the Lagrangian (H). They 
can be written as (recall that our metric tensor used to rise and lower spacetime 
indices is diag(l, —1), and that we use a summation convention for spacetime indices 
/i,z/G {0,1}) 

[D^,F'^'']=f, F^, = -[D„D,] (28) 

with = and j° = p. Here we introduced the covariant derivative 

= d,, + igA^. (29) 

In this paragraph (only!) we use a notation which takes into account the Leibniz 
rule for differentiation. In this notation one distinguishes, e.g., [d^A^) from dyA^ = 
{di,A^)+A^dy which is obtained by using the Leibniz rule. Thus {d^Afj) = [di,, A^]. In 
this notation, gauge transformations Eq. (^) can be simply written as U^^D^U 
and — >• U~^j^U, and it is thus obvious that the Lagrangian equations above are 
gauge invariant. These equations imply 

[Dy,f] = [Dy, [Z^^,Fn] = ^9[F,y,F^''] + [D„ [1^.,^^] = -[D„f] 

where we used the Jacobi identity, the anti-symmetry for the commutator, and 
[F^y^Ff""] = (cf. Eq. d)). We thus also have the equation [Dy,f] = [Dq, p] = 0. 

We can now forget that these differential equations originate from a Lagrangian 
and instead use them as definition of our Yang-Mills model. We find it convenient 
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to write these equations in the following form (we use the definition in Eq. (^; to 
simplify notation we write (df^X) as d^X in the sequel), 



doAi = E + diAo + ig[Ai, Ao] 



(30) 
(31) 
(32) 
(33) 



doE + ig[Ao,E] = 
dop + ig[AQ,p] = 
diE + ig[Au E] = p. 



Note that gauge invariance trivially implies the important 

Proposition: If An{t,x), E{t,x) andp{t,x) is a solution to the Eqs. j^S^j)-^^), then 
A^(t,x),E^(t,x) and p^{t,x) Eq. is also a solution for arbitrary differentiate 
GL^ -valued functions U(t,x) on spacetime. 

We now rewrite these equations in our two different gauges. 
Diagonal Coulomb Gauge 

As argued in Section ^ it is always possible to find a gauge transformation (given 
by Eq. (pH]) ) bringing Ai into the form Eq. (^. We now impose this gauge, i.e., we 
replace A^,E and p by A^,E^ and p^ with U given by Eq. (PT|). Using Fourier 
transformation, Eqs. (|30D -(p3|) can be written as (to simplify notation, we still use 
the symbols E, Aq, p instead of E'^ , Aq , p^) 



27r5„,o5"'''9og"(t) = ^"^(t,n)+^(n + (7[g°(t)-g^(t)])io^(t,n), (34) 

N 

doE-^it, n) = E E {^o\t, k)E^^{t, n-k)- E"\t, k)Af{t, n-k)), (35) 



9op°^(t, n) = E E {^o'{t, k)p^^{t, n-k)~ r{t, k)Af{t, n-k)), (36) 



and n = 0. We will use the following notation for the unspecified components. 



fcgZ 7=1 



N 



fcez 7=1 




Ar(t,0) :=a"(t) 



(38) 
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which are arbitrary functions. They correspond to the residual gauge freedom which 
is not fixed by the diagonal Coulomb gauge, as discussed after Eq. (JT^). 
Eq. {m yields 



Af{t, n) = a"(t)5"'^5„,o + (l - 5"^5„,o 
and Eq. (|37| ) implies 



in + g[q-it)-qP{t)]y 



(39) 



{n + g[q-{t)-qf'{t)]) 



(40) 



(we used the notation E°"^{t, 0) = 27rp"(t) introduced in Section |3.1|) . By combining 
these equations we obtain 



A, 



{t,n) = a"(t)5"^(5„,o + (l - S'^fSnfi) ^ 



p"^(t,n) 



n + g[q-{t)-qf'{t)]y 



(41) 



We now use these relations to derive the time evolution of our dynamical quantities 
in the diagonal Coulomb gauge. Putting P = a and n = in Eq. (^) and using 
again E'^"{t,0) = 27rp"(t) we obtain 



aog"(t)=p"(t). 



(42) 



The time evolution of the p°'{t) = E°"^(t,0)/2-K follows by inserting Eqs. ( ^OD and 
(IT]) in Eq. (H): 



2^ EE (i-^"Xo) ^"^(^'-)^^"(^'--) 



" (2vr)^ ^ (n + glq'^it) - g^(t)])^ 

Inserting Eq. (^) and (|4ll) in (|3^) we get the time evolution of p°'^{t,n) 



(43) 



n 



^p-^(t,n) (a"(t)-a^(t) 
Zn ^ 



(fc + ^[g"(t)-g7(t)])2 

X p'^'^{t,n-k)p'rP{t,k) 
^ '''Uk + 9[q'^it)-qm])' 



(44) 
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We now have obtained the time evolution equations for our dynamical system in 
a direct way. This result justifies our derivation of the Hamiltonian (p3D in the 
last section: The equations are identical with the Hamiltonian equations 

obtained from the Hamiltonian and the Poisson brackets (The proof 

of this is a straightforward computation which we skip.) 

It is worth to note that and are not necessarily real. Denoting the real and 
imaginary parts of by and pf, respectively, and similarly for g", we obtain from 
Eq. (H) 

5og"(t)=K(t), 



and from Eq. 



dop'^it) = Re 



29 
(2vr)2 



N 

EE( 

n ez 7=1 



p"^(t,n)p^°(t,-n) 



n,0) 



{n + g[q?{t)-qnt)+^[qnt)-qUtmr 



and similarly for the imaginary parts. Similarly, the time evolution of Eq. ( ^41) can 
be split in real and imaginary parts. In this general case, we can interpret the model 
as describing 2N particles interacting with each other and spin degrees of freedom. 



Weyl gauge 

We now discuss our dynamical system in the Weyl gauge Eq. (|T3|) . Imposing this 
condition, (pO|)-(|33|) are equivalent to 



doAi = E, doE = 0, dop = 
diE + tg[Ai,E]= p. 



(45) 



which can be solved trivially (see Section |^). We also observe that these equations 
are equivalent to the Hamiltonian equations obtained from the Hamiltonian ( ^2li) and 
the Poisson brackets 1 ^1^ and ^IBj). (Again the proof is a simple computation which 
we skip.) We thus have rigorously justified all results which in Section |3.1| using the 
Hamiltonian method. 



4 Solvable non-linear systems 

In this section we discuss how to obtain and solve dynamical systems which are special 
cases of the system defined by Eqs. (^)- (|4^) . We will impose additional constraints 
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on the charges p in order to obtain simpler systems which have natural physical 
interpretations and which can be solved explicitly. This will be done according to 
the following recipe: 

1. Make an ansatz for the charge x) which is parameterized by a finite number 
of X- independent functions and which is consistent with Eqs. (^Tj) and (Q). 



2. Insert this ansatz for p in the Eqs. (^21) -(^). To obtain the Hamiltonian and 
Poisson brackets giving rise to these equations, insert this ansatz for p in the 
Hamiltonian (^) and the Poisson bracket (p5|). 

We will also show how to obtain a one-parameter family of Lax pairs and conser- 
vation laws in our formalism. 

We are interested in the solution of the system of differential equations (^21) -(^) 
with the initial conditions, 

g"(0)=go", P"(0)=K, p°^(0,n)=pf(n). (46) 

Our method of solving this initial value problem is based on the results obtained in 
the last section. It was shown there that the Eqs. (^)-(|4^) are obtained from a Yang- 
Mills gauge theory by imposing the diagonal Coulomb gauge Eq. (^. Moreover, this 
very gauge theory in the Weyl gauge Eq. (|13D leads to the /mear differential equations 
(^Sp which can be solved trivially, 

E{t,x) = E{0,x), 
Ai{t,x) = Ai{0,x) + E{0,x)t, 
p(t,x) = p(0,x) (47) 

with the initial data £'(0,x), Ai(0,x) and p(0,x) satisfying Gauss' law Eq. (PBj). It is 
important to note that our solution (^) satisfies the Gauss' law for all t if it satisfies 
it for t = 0. From our discussion in the last section it is obvious that initial conditions 
for the Yang-Mills fields corresponding to the initial conditions Eq. (^) are, 

(0, x) = 6"^q^, r dxE"^{0, x) = 27ip^, p"^(0, n) = pf{n). (48) 



Gauss' law Eq. (|33D then completely determines E°'^{0,x) and thus the solution of 
our gauge theory model in the Weyl gauge. In Section ^72| , we also discussed how 
to construct the gauge transformation transforming to the diagonal Coulomb gauge. 
Combining these results we obtain the following recipe for solving a large class of 
non-linear integrable systems: 
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1. Use the initial conditions Eq. (pS]) to calculate E{0,x) from the Gauss' law Eq. 

2. Take Ai(t,x) given by Eq. (|4^) and calculate S(t,7i) from Eq. (||). 

3. The positions q°'{t) at time t are given by the eigenvalues A"(t) of S{t,iT) ac- 
cording to A"(t) = e-*^'^^''"^*), provided that A"(t') ^ A'^(t') for all a ^ /? and 
< t' < t. If this latter condition holds we say that no particle collisions occur 
in the time interval [0,t]. 

4. The solutions p{t,x) are given by U~^(t,x)p{0,x)U(t,x) where U{t,x) is the 
gauge transformation defined in Eq. (|Tl|). 



We note that the no-collision condition above is due to the regularity assumption 
which we needed to prove the existence of the diagonal Coulomb gauge in Section 



2.2. 



We now show how to obtain Lax equations and conservation laws in our for- 
malism. For that we observe that one of our Yang-Mills equations, d()E{t, x) + 
ig[AQ{t,x), E{t,x)] = 0, has the form of a Lax equation for arbitrary fixed x: For 
any x, L{t) = E{t,x) and M{t) = igAQ{t,x) is a Lax pair, dL/dt + [M,L] = 0. We 
thus have a one-parameter family of Lax pairs. By a standard argument this implies 
that tr[E{t, x)"] is invariant under time evolution, for arbitrary positive integer n and 
X G [— 7r,7r]. We have already found these conservation laws by a different argument 
at the end of Section 

In case the computation here can be done explicitly we obtain explicit solutions. 
Lax pairs and conservation laws. We will now discuss several examples of increasing 
complexity where this is possible. Before that we state three identities which will be 
important for us. 



Three identities. The first identity we will need is well-known, 

El vr^ .X 

7 = — 9-7 — T, r eC (49) 



see e.g. fl^). The second is a generalization of the first 



^7 — ^x = e '"^M— ^— — + z7rs2^cot(7rr)-7r|s2^| , r G C, s G R, (50) 
^ (n + r)2 \^sm (vrr) / 
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where 

S2n '■= s — 27m with n G Z such that — tt < s — 27m < n. (51) 

Finally, 

E — = ^^-^^2. sGM. 52 

;^^o 2 3 

Proofs of Eqs. (^Dj) and (|52|) are given in Appendix ^ 



4.1 Example 1: Spin CM model with v{r] 



,-2 



As the simplest example we now show how to obtain and solve the Calogero model 
and certain spin generalizations thereof. 
For that we make the following ansatz, 

p°^(t,n) = 5„,o(l - 5°^)27r(7s-^(t). (53) 

It is important to note that Eq. (|^) implies that dQp°'^{t,n) = for n ^ and 
(?op""(t, 0) = 0. Hence the ansatz Eq. (|53|) is consistent, and the functions s"''^(t) for 
a ^ f3, together with the p"(t) and q°'{t), are the dynamical variables of a dynamical 
system. 

From the Hamiltonian (^) and the Poisson brackets (|2^ ) -(|25| ) we obtain 

1 1 „o>l3„/3a 

= _ s"'/55/3'°) (56) 

which determines the equations of motion for the dynamical variables. Moreover, the 
initial conditions are 

g°(0) = q^, p"(0) = and s"^(0) = sf all real. (57) 

We now have to determine the corresponding initial conditions for the Yang-Mills 
fields. Using E"^(0,x) = ^ E„ez ^°^(0, n)e*"^ and Eqs. (^) and (|Dp we obtain 



i?"^(0, x) = p-S-f' + (1 - r^) ,J^° (5^ 



a/3 
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Moreover, we recall Af{0 ,x) = Inserting this and Eq. (|58| ) in Eq. ( |4^ ) we 

find that Ai(t,x) = Ai(t) is independent of x, 



Af{t) = q^6^^ + ( + (1 - 5"^)—^ 1 t. (59) 



q^d""^ + + (1 - 5°^ 




This is the solution of the Yang-Mills equations. To find the solution for our dy- 
namical system we need to transform this to the diagonal Coulomb gauge Eq. (|^). 
This is simple here: all we need to do is diagonalize the matrix Ai[t). Furthermore, 
our discussion in Section ^ implies that s°'^(t) can be obtained from the matrix U(t) 
diagonalizing Ai{t). We thus obtain the following 

Result 1: Provided that no particle collisions occur in the time interval [0,t], the 
solution q°'{t) of the initial value problem corresponding to the Hamiltonian and 
the Poisson brackets (|5^ and ( |5^ j is given by the eigenvalues of the matrix Ai{t) in 



Eq. (Hi; 

f/(t)-%(t)[/(t) = diag (g^(t), g2(t), . . . , 

Moreover, the diagonalizing matrix U {t) gives the time evolution of the spin degrees 
of freedom, s{t) = U{ty\s{0)U{t). 

We thus have recovered the integrable systems and their solution previously found 
in Ref. In the present case we only get one Lax pair (the Lax pair family is 
independent of x) which coincides with the one given in Ref. Q]. 

Note that the special case where all initial s"^ are identical, s"^(0) = e, corre- 
sponds to the Calogero model, i.e., the Hamiltonian Eq. (^ with v{r) = r~^. To see 
this, it is important to notice that the functions a"(t) in Eq. (|^) can be chosen such 
that s"^(t) = e for all t. Indeed, with 

27re ^ 1 
a-(t) = -— y , ,,,, (60) 

9 (g"(t)-g^(t))' ^ ' 

we obtain 5os"^(t) = 0. 

Remark: It would be interesting to find and study a case where collisions occur. 
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4.2 Example 2: CM model with v{r) = a^sin ^(ar) 

In this section we show in detail how to recover the known solution of the Sutherland 
model, i.e., the Hamiltonian Eq. (|l]) with the potential v{r) = a^sin~^(ar). This is a 
simple special case and warm-up for what is discussed in Section 

To obtain the Hamiltonian of the Sutherland model we make the ansatz 

^"/^(O, n) = (1 - r^)27r^e Vn G Z (61) 

equivalent to 

p°/3^0, x) = (1 - S'"^)27igeS{x). (62) 

Choosing 

"°<''^?.&'-*°^*-'(™y^ (63) 

implies dop°''^(t,n) = 0, as can be seen by inserting Eq. (^) in Eq. (PI). We thus 
obtain p"^{t, n) = {l — 6°'^)2'n'ge for all t, and and are in fact the only dynamical 
variables in the system. 

Inserting Eq. ( |6lD into the Hamiltonian Eq. (123|) and using Eq. (^) we obtain 

2 ^^^sm (7r^[g"-g/^]) 

which for g = a/n equals the Hamiltonian of the Sutherland model as discussed in 
Section |l|. 

We now show how to solve the initial value problem for the Hamiltonian equations 
following from the Hamiltonian Eq. (|6^) and the Poisson brackets {q'^,p^} = S"^. We 
first have to determine the initial conditions for the Yang-Mills fields corresponding 
to g"(0) = Qq and p"(0) = all real. For t = we can write Gauss' law Eq. (|33|) as 

di (e^f^«^E(0, x)e-^^^''^) = e^^^°^p(0, x)e-^^'3ox ^g^^ 

where Qq = Q{0) = diag(gQ, . . . , g^) and p(0, x) is given by Eq. (|6^). Since p(0, x) = 
except for x = 0, it follows that e^^^°^E{0,x)e~'^^^'^^ is constant in the intervals 
— TT < X < and < x < tt. We therefore can write 

E{0, x) = e-^»««^5±e^^^«^ for x < (66) 
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with B± constant matrices. 

To determine the matrices B± we integrate Eq. (|65| ) from — e to + e and then 
take the hmit e | 0. This gives 

Bf - 5°^ = (1 - 5"^)27r^e. 

We also recall J^^dxE°"^{0,x) = 2ttpq which yields 

Moreover, the boundary condition £'(0, — tt) = E{0, ir) implies e^*'^^'3oi?_e~^*'^3'3o = 
-B+, or equivalently 

A straightforward computation shows that these equations determine the elements 
of the matrices -B_ and -B+ as follows, 

Bf = + 1 - ; (67) 

Inserting the initial condition Ai{0, x) = Qq and Eq. ( |66D in Eq. ( ^Tf ) we finally obtain 
the explicit solution of the Yang-Mills equation in the Weyl-gauge, 



> 



Ai{t, x) = e-^^^«^' {Qo + B±t) e^^^»^ for x < (68) 
(we used Qo = e-'^QoxQ^^igQoxy 

We now determine the gauge transformation which transforms this solution to the 
diagonal Coulomb gauge. For that we need to solve Eq. (||) for S{t,x). We define 
S{t,x) = e'3QoxS{t,x) and observe that Eqs. (|) and (pi) imply 



diS{t, x) + igB±tS{t, x) = ior x < 0. (69) 
This latter equation can be solved trivially, 

S{t, x) = e-*3^±*^'5(t, 0) for a; < 0. (70) 
We thus obtain S{t, vr) = e-'^'^^^+^Sit, 0) and S{t, 0) = e-^'^s^-*^(t, -vr), and 
S{t, tt) = e-'''^^''S{t, tt) = e-^"5'3og-*-<;B+t^(i^ g) = 

= ^-^'^dQo^-iTrgB+t^-ingB-t^-ingQo 
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where we used S(t, — vr) = /. 

In the present case we can easily prove that the Yang-Mills field Ai{t, x) Eq. (|6^) 
is generic, i.e., if Qq ^ for all a 7^ /3 then particle collisions never occur. This 
follows from the fact that the Hamiltonian Ti in Eq. (|6^) is conserved under the 
time evolution, and this implies an upper bound on sin~^(7r(yf[g° — q^])- There exists 
therefore an e > such that \q'^{t) — q^{t) — n/g\ > e for all a 7^ /5, t G M, and 
all integers n. Our discussion in Section |2.2| implies that the latter is equivalent to 
Aiit^x) being generic, i.e., regular at all times t. 

From the discussion in Section ^ we thus obtain the 

Result 2: The solution of the initial value problem for the Sutherland model Eq. 
(I^j is given by the eigenvalues A"(t) of the matrix 

with the matrix elements of B± given in Eq. and Qq = diag(gQ, . . . , q^), accord- 
ing to0 q^{t) = -\og{\^{t))/{2mg). 

Remark: To compare this with the well-known result reviewed in ^ we observe that 
E± = e±*'^»<5o5eTi7rgQo^ ^here 



isin(7r^[g^ - q^]) 

Using this it is easy to see that the matrix in Eq. (|7T1) has the same eigenvalues as 
the matrix 

^-ingQo^-2ingBt^-iTTgQo ^^-^2) 

which shows that the result above is identical with the one derived in Ref. 0. 

We finally show how to recover the Lax pair and conservation laws for the Suther- 
land model in our formalism. We observe that our computation of E{t, x) at t = in 
terms of = = 0) and = p°'{t = 0) immediately generalizes to t 7^ 0. From 
our solution at t = we therefore can read off that 

E{t, x) = e-*9{^'T-)«W5(t)e^f(^^")<3(*) for x < (73) 

with 



'*the branch of the log is fixed by continuity of q" (t) 
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Inserting Eqs. (0) and ( p^ ) in Eq. (^) and using Aoit, 



x] 



-y 



2n /-^nez 



Mo(t,n) to- 



gether with the identities in Eqs. (|50| ) and (|5 

/ 



eg 



AT 



E- 



TT 



we obtain after a simple computation 

A 



It sm\7Tg[q'^{t) - q^t)]) 



X 



TT 



^sin2(7r(7[g"(t) 



mxa. cot(7r(7[g"(t) -g^(t)]) 



We thus have found a one-parameter family of Lax pairs for the Sutherland model. 
However, due to cyclicity of trace, the corresponding conservation laws are indepen- 
dent of x: tr[E{t,x)^] = tr[B{t)"']. It is interesting to note that the standard Lax 
pair for the Sutherland model, derived e.g. in Ref. 0, is recovered from our Lax pair 
family above by setting x = n. 



4.3 Example 3: CM model with v{r) = a^sinh ^(ar) 

We now make a simple but important observation: In Example 2 above we can allow 
for general complex initial data and and thus obtain time evolution equations for 
complex valued functions q°'{t) and p°'{t). In this case all our discussion in Example 
2 above goes through, except the argument showing that no particle collisions occur. 
Thus our Result 2 above holds also for the general complex case at least locally in 
time, i.e., until the first particle collision occurs. 

An important special case is obtained if e is real and all q^ and are chosen 
purely imaginary, i.e., they are real if multiplied with i. It is easy to convince oneself 
that the results for this case we can obtain from the results in Example 2 simply by 
the following replacements 

-> p»{t)^zp^{t) (75) 

where and p°'{t) real. The Hamiltonian and Poisson brackets then get minus 

signs which can be conveniently removed by the transformation, 

n^-n, l,-}^ -{■,■}, (76) 
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which leaves the Hamilton equations doX = {Ti, X} invariant. The resulting initial 
value problem thus is equivalent to the one obtained from the Hamiltonian 



n = ^ Up"? + t E ..^/fi (77) 

2 ^1 2 smh {7Tg[q» - ql^]) 

and the standard Poisson brackets = 5"^. This initial value problem is iden- 

tical with the one for the Sutherland model with the potential v{r) = a^sinh~^(ar) 
and a = vrgf 0. In this case we can argue as in Example 2 that no particle collisions 
can occur, i.e., |g"(t) — 5'^(^)| > e for some e > and allt G M and all [3. Result 
2 above then implies: 

Result 3: The solution of the initial value problem for the Sutherland model defined 
by the Hamiltonian in Eq. ( fTTj ) is given by the eigenvalues A"(t) of the following 
matrix, 

with 

enge^ 



smh{TT g[q^ - q^\) 
and Qq = diag(gQ, . . . ,qQ), according to q^{t) = \og{\°'(t)) / {2TTg) . 

As in Example 2 one obtain Lax pair etc. and can check that what we obtain is 
equivalent to the known result for this model given in M. 



4.4 Example 4: Generalizing the Sutherland models 

We now discuss a general case for a spin-particle dynamical system equivalent to a 
Yang-Mills gauge theory as discussed in Section |] and which we can solve explicitly. 
Our result here generalizes the one in Ref. 

We choose the charge distribution of the form 

m 

p'^^{t,x)^Y.pfm^-^j) (80) 

where m is some positive integer and we assume, without loss of generality, 

Xo = — TT < Xi < X2 < . . . < Xm < Xm+1 = VT. (81) 
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This charge distribution describes matter charges locahzed at the points Xj, 1 < j < 
m. Using Eq. (p2D it is easy to see that this ansatz is preserved under time evolution. 
Moreover, using 

m 

p-^t,n) = Y.pf(t)e-^''''' (82) 

k=l 

we see that Eq. (^) imphes the constraint 

m 

EPr(0) = Va = l,...,iV (83) 

on possible initial conditions for the pj^. 

Using the identities ([50| ) and (0), the Hamiltonian in Eq. (|23|) becomes (we 
specialize to a'^(t) = for simplicity) 

-\ N -i N m -I N m 

n = - E(p")' + 2 5: E ^A<i'' - <f)pfpT + 2 ^ ^ ^^prpf' (84) 

a=\ a,0=lj,k=l a=lj,k=l 



where 



and 



vAr) = \e-^'"^^^ ( + ^ cot(vr^r) - ^ ) (85) 



(86) 

^ S7r2 47r 12 ^ ' 



Here Xjk is defined as 

= {xj — Xk)2TT- (87) 
The Poisson brackets for this system are: 



{pf , = ^^2vr5,, {S^-'pf - S^'yf) . (89) 

We now solve the equations of motion for the system given by the Hamiltonian 
(|84D and the Poisson brackets (|88|)-(|89D. More general than in [|^, we do not assume 
that q°'{t) etc. all are real. 

We start with Gauss' law Eq. ( p5| ) where Qq = Q{0) is given by (|^ and p{0,x) 
by (|80|) . Since p(0,x) = except for x = Xj, we obtain from that E{0,x) = 
Q-igQox j^^igQox ^ when Xj < X < j = 0, . . . , m, where Bj is some constant matrix 
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(we recall that xq = — vr and Xm+i = tt). To determine the matrices Bj, we integrate 
Eq. (|65D from Xj — e to Xj + e and then take the limit e | 0. This gives the recursion 
relations Bj - = e'^^^''^ pj{0)e-'3Qoxj ^ ^^^^ ^-^e condition E(0, -tt) = E{0,tt) 

implies e^*^'^o'^i?oe~^*^'^"'^ = B^- Combining these equations and using the definition 
Eq. (|4^ ) of p", we obtain after a straightforward calculation, 

p"(o) + Epr(o) - E ^^^^^ Epr(o) + 
^=1 i=i ^=1 / 

^ V 2zsin(^?7r[g"(0)-g/^(0)]) ^ ^ 

with sgn(x) = 1 for x > and —1 for x < 0. With that, the solution of the Yang-Mills 
equations in the Weyl gauge, Eq. (|47|), gives 



Ai(t, x) = e-*3'3ox ^ ^^.^^ gigQo^ for Xj<x< Xj+i. (91) 

Next, we solve diS{t,x) + igAi(t,x)S(t,x) = which is equivalent to 

diS{t, x) + igBjtS(t, x) = for xj < x < xj+i (92) 

for S{t,x) = e^3Qo^^(t,x). This implies S{t,x) = e-^3^^*(^-^^)^(t, x^) for Xj < x < 
Xj+i, thus 

S{t,Xj+i) = e-*9Qo-^+ie-*3^^*("^+^-"^)^(t,x,-) = . . . = (93) 

_ g^-'tgQoXj+i^-igBjt{xj+i-Xj)^-igBj^it{xj-Xj-i) _ _ _ ^-igBot{xi+Tr) ^-igQoTr 



where we used S{t, — vr) = /. Especially (for j = m), 

_ ^-igQoTT Q~igBmt{n-Xm) . . . 
y.^-igB^^it{xm-XrTT,-i) _ _ _ ^-igBot{xi+n) ^-igQoTT ^g^-j 

From our discussion in Section ^we can obtain the q"{t) from the eigenvalues of this 
matrix. Moreover, 

p,(t) = f/(t,x,)-V,(0)f/(t,x,), 
where U{t,Xj) is given by (|TT]). We thus obtain the 



Result 4: Provided no particle collisions occur in the time interval [0,t], the solution 
q"'{t) of the initial value problem defined by the Hamiltonian (j^[) and the Poisson 
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brackets is given by the eigenvalues A"(t) of the matrix S{t,TT) defined in 

Eqs. and ^9di), according te^ q°'(t) = —\og{X"'{t))/{2Tcig). Moreover, the time 

evolution of the spin degrees of freedom is given by 

p.{t) = e-*3(^'^+-)«W\/(t)-i5(t,Xj-)"Vj(0)5(t,Xj-)^We*^^^''^"^^^*^ (95) 
where S(t, xj) is defined in Eqs. ( ^ ) and ^^), Q(t) = diag(gi(t), q2(t), . . . qN(t)) and 



-ig2-KQ{t) 



V{t) is the matrix diagonalizing S{t,TT), i.e. 

V{t)-^S{t,7T)V{t) 



As discussed in the beginning of this section, tr[E{t,x)^] is conserved under time 
evolution. For the special choice of charges studied above, we can evaluate E(t,x) 
and obtain E{t,x) = e~'^^'^^*^^ Bj{t)e^^^^^^^ for Xj < x < Xj^i where Q{t) and Bj{t) are 
as in Eqs. (0) and <^ but with g"(0), p"(0), and p"'^(0) replaced by p"(t), and 
pj^{t), i.e., the solution of the initial value problem which we solved above. Using 
cyclicity of the trace, we conclude that tr[Bj{t)"'], for an arbitrary positive integer n 
and j = 1, . . .m, are time independent: Each of them is a conservation law. 

The one-parameter family of Lax equations is doE(t, x) +ig[AQ(t, x), E(t, x)] = 0. 
We now compute the corresponding family of Lax pairs explicitly, similarly as in 
Section |4.2|. We obtain after a straightforward calculation. 



(J 7T1, t \ 

e=i i=i £=1 / 

V J 2tsm{gn[q-it)-q^it)]) ^ ' 



and 



Af{t^x) = -5-^^ tpTit) - vr|(x - x,U + ^) + 

2,tl ' ^ ^ \sm^ing[q'-it)-qPit)]) 

i{x — Xk)2TT ./ r a/,\ I3u\}\ I ~ ^A:)27r| 1 



+ '-^^^—^cot{ng[q-{t)-q^m 

TT 71 



^the branch of the log is fixed by continuity of q" (t) 
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Note that for m = 1, the dynamics of the spin and the particles decouple, and we 
recover the results in Sections ^]2| and [4.3| . 

Note that all variables q°'{t), p°'{t), and p^^if) above were allowed to be complex. 
If we restrict them to be real we recover the results in 0]. This generalizes the results 



for the Sutherland model is Section [4.2| . If we assume that all q°'{t) and are 
purely imaginary and the pj^it) all real we obtain results generalizing the ones in 
Section ^j3| (similarly as discussed in Section ^Tsj they can be obtained from the ones 
for the real case by the substitutions in Eqs. ([75|) and ([761). 

4.5 Example 5: Novel systems 

The previous examples all where such that it was possible to give very explicit solu- 
tions. As a final example we present a case where the solution can be only described 
in a somewhat implicit way. It corresponds to novel integrable spin-particle systems 



which generalizes the model discussed in Section 4.1. The arguments here are similar 



to the ones in Section 4.4 and therefore we are rather brief. 



p"^(t, x) = ^e-'^[''"W-''''W]"sf (t), Xj.i <x<Xj (97) 



We chose the charge distribution as follows 

1 

where Aj = {xj — Xj-i), j = 1, . . . ,m, and 

Xo = — vr < xi < X2 < . . . < Xm-i < Xm = t: (98) 

with m some positive integer. This ansatz describes (essentially) piecewise constant 
matter sources and is preserved under time evolution (cf. Eq. (^)). Moreover, using 

p»^[t,n) = J2sf{t)— I da;e-*^[^"W-^^W]"e-'"" (99) 

j=l "'^'i-i 

we see that Eq. (^) implies the following constraint 

m 

^sj"(0) = Va = l, ...,iV (100) 
i=i 

on possible initial conditions for the s"^. 
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Using the identities ( ^0] ) and (0), the Hamiltonian (pSf ) becomes (again we set 
a"(t) = for simphcity) 

1 N T N m -I N m 

n = 2 Eb")' + 2 ^ E ^A<i'' - /)^f + 2 ^ ^ ^^-^^^"^r (101) 

a=l Q,/3=l7,fc=l «=lj,fc=l 

where 

X f J , + ^^"-^^^- cot(.,r)- '("-^^^-' ') (102) 

\^Sin (TTf/r) TT 71 J 

and 



(of course these integrals could be further evaluated). The Poisson brackets for this 
system are: 

{sf,si''} =^g27r6,,{6^'^'sf-6^''^sf^) (105) 
where the last equation is obtained by inserting 

sf{t)= r cixe^3['?"W-^''Wl^p"^(t,x) (106) 

in Eq. (|18D. 

We now show how to solve the initial value problem for the Hamiltonian equations 
of this model. Writing Gauss' law in the form Eq. (|65D and introducing (Qo = <5(0), 
as in the previous Section) B(0,x) = e^^'^°^ E{0, x)e~^^^°^ we obtain with Eq. (P7D, 

E"/^(0, x) = Cf{0) + j^^-^^sf (0) for <x<x, (107) 

(Xj-Xj_i) ' 

where 

i-i 

Cf (0) = 5°^(0, -vr) + 5] sf (0). (108) 

k=l 
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Similarly as in Section |4l^ , the conditions 27rp"(0) = J^^dx E°"^{0,x) and E{0,it) 
E{0, — tt) determine B{0, — vr), and we obtain after some computations, 

/ ™ A- / 

5-/^(0, -tt) = rMj9"(0)-E^UsJ"(0) 

+ (l-*")2?ii«oF^g4^(0). (109) 

We thus obtain Ai{t, x) = e'^^^"^ [Qq + 5(0, x)] e'^Qo^. 

To find the gauge transformation bringing this to the diagonal Coulomb gauge we 
define S{t, x) = e^^^'^^S{t, x) converting the differential equation in (||) to diS{t, x) + 
igB{0, x)tS{t, x). Requiring S{t, —n) = / we obtain 

S{t,x) = e-^Q°%- (t, 1^^) T,-i(t, 1) ■ ■ ■Ti(t, l)e~^sQon < x < x, 

(110) 

where 

T,(t,r) =Vexp (^-igtA, J^^ d^[Cj{0) + ^Sj{0)]^ . (Ill) 

This gives the solution of the model: The q°'{t) are determined by the eigenvalues 
of S{t,7r) as in the previous examples. To compute Sj{t) is more complicated here, 
however: One can use Eq. ( |106| ) with p{t,x) = U^^{t,x)p{0,x)U{t,x) and U{t,x) 
determined by S(t,x) as discussed in Section ^]2[ In the present case we could not 
make the solution more explicit. 

It is possible to write down explicitly the Lax pair family for the present model, 
and there is also a 1-parameter family of conservation laws here, namely tr[B{t,x)"'], 
n e N and x G [-vr,7r], where B{t,x) is as in Eqs. ([107|)-([109D with g"(0), p°(0) and 
Sj(0) replaced by g"(t), p°'(t) and Sj(t). 



A Appendix: Proofs 

In this Appendix we prove the identities Eqs. (^) and (^2|). To prove Eq. ( pOD we 
define the following function. 



pi(n+r)s 
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(r is a complex and s a variable). We observe 

Integrating this twice gives 

g{r, s) = A{r) + B{r)s — vr|s| for — vr < s < vr (113) 

where A{r) and B{r) are integration constants. We now determine A[r) and B{r). 
Using Eq. (|49|) we conclude that 



1 _ TT^ 

('^ + ''^y sin^vrr' 



Moreover, the definition Eq. (|1 1 2|) of g{r, s) implies 



Inserting Eq. (|113|) we get 

B{r) = in cot (vrr). 

We thus obtain, 



V -r = e — 2 ^ ^^-5 cot (yrr) - 7r|s| for -vr < s < vr. (114) 

'n + r)"* \sin vrr / 



The l.h.s. of this equation obviously is periodic in s with period 27r. To extend the 
r.h.s. to all real s we therefore only need to replace s by := s — 2?™ with the 
integer n chosen such that — tt < s — 27m < n. This proves Eq. (0). 
To derive Eq. (^) we define 



gins 



(s is real) and compute 

^i-^ = _ y e'ns ^ _2TrS(s) + 1 for -vr < s < vr. 
ds t^z 
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Integrating this twice yields 



S2 



h{s) = -n\s\ + — + Cs + D 



where C and D are integration constants. Periodicity of h{s) imphes C = 0, and the 
formula 

°° 1 _ 

6 



n=l 

2 



see, e.g., |jT2|) determines D = Thus 



— Y = + ~ + ~ ~^ < s < TT. (116) 
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Again the r.h.s. can be extended to all real s by replacing s by S2n as defined above. 
This proves Eq. (|52D. 

B Appendix: Poisson bracket for the charges 

In this appendix we derive the Poisson bracket for the charges p{x), i.e., Eq. ([T8|) . 
As mentioned in the text, consistency of Gauss' law G{x) ~ requires 

doG"'^'{y) = {G"'^'{y),H}c^O Va', /?' = 1, . . . , iV. (117) 

Inserting Eqs. ([ISD and ([17|) we can rewrite {G"''^' (y), H} ~ as 

N 

- / dx 

27 



^ dx Y: A^,-{x){p-^{x),p-'^'( 

1 /-TT ^ 

a,/3,7=l 

- ^^"(x) {a.E°'^(x),<^(2/)E^^'(y) -E"'^(y)Af' (I/)} 

- A^,"{x) {Ar{x)E^^{x) - E"^(x)Af (x), 9,E"'^'(y)} - tgA^o^i^) 

X {Ar(x)i?^^(x) - E'^\x)Af{x),A^''^{y)E^P\y) - E'^'^{y)Af (y)}) . 
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Here we used that the only non-vanishing Poisson brackets are those for Ai with E 
and for p with itself. Using now the Poisson bracket (|T6D the r.h.s. of this becomes 

1 rn N 

- dx Y: A^,-{x){p-^{x),p-'^'{y)} = 

•^-'^ a,/3=l 

= + tg dxY (A^'\x)E^^'iy)dJix - y) - Af {x)E'''\y)dJ{x - y) 

-'"^ 7=1 

+ Af{x)E^'^{x)dJ{x - - A''^\x)E^P' {x)dj{x - y)) 

- {tgf £ dx 5{x - y) {[[A,{y). ^0(2/)] , E{yr'^' + [[Ao(y), E{y)] , AMT'^') 
= zg dx 5{x - y) (- [d^E{y), Ao{x)r'^' - ig \\A,{y), E{y)\ , Ao{x)r'^') 

= i-ig) r dx6{x-y)[^^E{y)+^g[A^{y),E{y)],Ao{x)f''' 

J —IT 

- i-ig) f dx6{x -y)[p{y),Ao{y)f''^' 

J —IT 

= ± dx J2 

Here the second equality follows by using partial integration and the Jacobi identity, 
and in the fourth equality Gauss' law was used. From this the Poisson brackets in 
Eq. (|18D follow. 
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